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This  study  utilizes  the  unique  data  set  obtained  during  the  Structure  and  Atmospheric  Turbulence 
Environment  (STATE!  experiment,  conducted  during  June  1983  at  Poker  Flat.  Alaska,  to  examine  the 
structure  and  characteristics  of  the  wave  field  near  the  summer  mesopause.  It  is  shown  that  the  rocket 
and  radar  data  together  permit  a  much  more  detailed  specification  of  wave  parameters  than  would  be 
possible  using  either  data  set  alone.  The  results  of  this  analysis  suggest  that  the  wave  field  near  the 
summer  mesopause  is  composed,  in  general,  of  a  superposition  of  wave  motions  which  act  collectively  to 
produce  regions  in  which  the  wave  field  is  convectively  or  dynamically  unstable.  These  regions  are  found 
to  correlate  well  with  zones  of  enhanced  turbulence  and  small-scale  wave  activity,  suggesting  the  pro¬ 
cesses  and  elVects  of  wave  field  saturation. 


1.  Introduction 

There  has  been  a  resurgence  of  interest  in  the  last  few  years 
in  the  propagation  and  saturation  of  gravity  waves  in  the 
middle  atmosphere,  following  the  recognition  of  their  role  in 
establishing  the  large-scale  circulation  and  thermal  structure 
of  this  region  {^Houtfhion.  1978;  Lindzen.  1981].  Initial  mod¬ 
eling  studies  by  Matsuno  [1982],  Holton  [1982.  1983],  and 
Dunkerton  [1982a]  showed  that  crude  parameterizations  of 
gravity  wave  drag  and  induced  diffusion  could  account  for 
several  significant  features  of  the  middle  atmosphere  circu¬ 
lation  and  structure.  These  include  the  reversals  of  the  vertical 
shear  of  the  mean  zonal  wind  and  of  the  mean  meridional 
temperature  gradient  in  the  mesosphere.  Similar  results  \.v.re 
obtained  in  a  high-resolution  study  by  Mi\ahara  el  il.  [1986] 
that  employed  direct  simulation  of  gravity  wave  motions. 
Other  theoretical  studies  addressed  several  ispects  of  gravity 
wave  propagation  and  saturation  of  relevance  to  »his  observa¬ 
tional  study,  including  effects  of  wave-wave  and  ave-mean 
flow  interactions  [Wein.slock,  1982;  Dunkerton,  1982/>;  Wal- 
terscheid,  1984;  Fritts,  1985],  convective  adjustment  of  unsta¬ 
ble  wave  fields  [Dunkerton  and  Fritts,  1984],  and  the  induced 
diffusion  of  heat  and  constituents  {_Schoel  jrl  et  ai,  19,23;  Chao 
and  Schoeherl,  1984;  Fritts  and  Dunkerton.  1'’85;  Strohel  el  a!., 
1985], 

In  contrast  to  the  relatively  large  number  of  theoretical 
studies  of  gravity  wave  propagation  and  saturation,  only  a  few 
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observational  studies  have  been  performed  that  specifically 
address  gravity  wave  saturation  and  its  effects.  A  measure  of 
the  zonal  drag  due  to  saturating  gravity  waves  as  well  as 
phase  speed,  period,  and  wavelength  estimates  were  provided 
by  yinceni  and  Reid  [1983],  Related  studies  by  Smith  and 
Fritts  [1983],  Baisley  et  al.  [1983],  Meek  et  at.  [1985a],  Smith 
et  al.  [1987],  Frills  and  Vincent  [1987],  and  Reid  and  Vincent 
[19870.  h]  provided  information  on  gravity  wave  scales,  am¬ 
plitudes.  and  fluxes,  as  well  as  the  conditions  leading  to  satu¬ 
ration  of  the  gravity  wave  spectrum.  Additional  evidence  of 
saturation  and  of  characteristic  wave  amplitudes  and  their 
variability  were  presented  by  Vincent  [1984],  Meek  et  al. 
[19856],  and  Vincent  and  Frills  [1987J.  The  initial  theoretical 
and  ob.servational  studies  of  gravity  waves  in  the  middle  at¬ 
mosphere  were  reviewed  by  Fritts  [1984a]. 

Most  previous  studies  of  gravity  waves  in  the  middle  atmo¬ 
sphere  have  used  data  from  a  single  observing  system.  This 
usually  provides  inlormalion  on  only  a  single  variable  (veloci¬ 
ty,  temperature,  or  density)  and  causes  some  ambiguity  in  the 
description  of  the  gravity  wave  field.  When  two  or  more  in¬ 
strument  systems  are  used  together,  however,  a  much  more 
complete  determination  of  the  wave  field  and  its  middle  atmo¬ 
sphere  effects  can  be  obtained.  Such  was  the  case  in  the  Struc¬ 
ture  and  Atmospheric  Turbulence  Environment  (STATE)  ex¬ 
periment.  which  provides  the  basis  for  this  series  of  papers. 

The  simultaneous  and  nearly  colocated  radar  and  rocket 
data  collected  during  the  STATE  experiment  provide  a  unique 
opportunity  for  a  detailed  look  at  the  role  played  by  gr;  vity 
waves  in  middle  atmosphere  dynamics.  As  will  be  seen  later, 
this  data  facilitates  the  inference  of  gravity  wave  scales,  direc¬ 
tions  of  prop  gation,  mechanisms,  and  likely  consequences  of 
dissipation  as  well  as  the  middle  atmosphere  response.  Of 
particular  interest  in  the  STATE  experiment  are  an  expla¬ 
nation  for  the  very  strong  MST  radar  echoes  observed  near 
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POKER  FLAT  MST  RADAR 
15  June  1983 

HORIZONTAL  WIND  VECTORS 

(1 5-Minute  Averaged  Value) 


Fig.  I.  Time-height  cross  section  of  the  inferred  horizontal  wind  field  for  June  15.  1983.  Each  vector  represents  a  15-min 
average.  Note  the  low-frequency  motion  dominating  the  last  half  of  this  interval. 


the  summer  mesopause  and  the  mechanism  or  mechanisms 
that  provide  a  large  energy  dissipation  rate  (k)  and  a  small 
Kolmogorov  inner  scale  (/„)  in  that  region.  Estimates  of  these 
quantities  during  the  STATE  experiment  are  provided  by  Wat¬ 
kins  et  al.  [this  issue],  Kelley  et  al.  [this  issue],  and  Ulwick  et 
al.  [this  issue].  In  addition,  this  data  will  help  us  to  assess  the 
validity  of  the  various  theories  that  have  been  advanced  to 
account  for  previous  gravity  wave  observations  and  their 
middle  atmosphere  effects. 

An  overview  of  the  STATE  data  used  in  the  wave  field 
analysis  is  presented  in  section  2.  A  significant  feature  of  these 
data  is  the  very  high  degree  of  spatial  and  temporal  varia¬ 
bility.  An  elementary  discussion  of  gravity  wave  structure  and 
of  the  relationship  between  the  various  wave  parameters  ob¬ 
served  with  the  radar  and  rocket  instrumentation  is  provided 
in  section  3.  This  suggests  some  redundancy  in  our  ability  to 
estimate  various  gravity  wave  parameters.  The  analysis  of  the 
wave  field  observed  durirp  salvo  7  js  presented  in  section  4. 
These  Oatu  provide  clear  evidence  of  both  large-  and  small- 
scale  wave  motions  throughout  the  height  interval  over  which 
data  were  collected.  The  large-scale  motion  is  observed  in 
both  velocity  and  temperature  fields  and  is  seen  to  correlate 
well  with  the  radar  signal-to-noise  (S/N)  data  and  regions  of 
enhanced  small-scale  wave  amplitudes,  suggesting  that  insta¬ 


bility  and  dissipation  are  highly  localized  within  the  wave 
field.  A  second  event  yielding  a  similar  correlation  is  described 
in  section  5.  The  conclusions  of  this  study  are  presented  in 
section  6. 


2.  Data  Overview 

The  objectives  of  the  STATE  experiment  and  discussions  of 
the  rocket  and  radar  data  collected  are  presented  in  other 
papers  in  this  series.  The  purpose  of  this  section  is  to  review 
that  portion  of  the  data  that  can  be  used  to  infer  the  charac¬ 
teristics  and  effects  of  observed  internal  gravity  wave  motions. 

Perhaps  the  most  striking  feature  of  the  radar  data  collected 
near  the  summer  mesopause  at  Poker  Flat  is  the  enormous 
temporal  variability  of  the  velocity  field.  The  time-height  cross 
sections  of  the  15-min  average  inferred  horizontal  wind  be¬ 
tween  80  and  90  km  are  shown  for  June  15  and  17,  1983,  in 
Figures  1  and  2.  The  first  of  these  exhibits  relatively  high- 
frequency  motions  (T  S  2  hours)  at  upper  levels  from  0000  to 
0800  LT,  Thereafter  the  motion  field  is  dominated  by  a  low- 
frcquency  wave  (T  ^  7  hours),  which  rotates  anticyclonically 
with  time  and  height,  suggestive  of  the  low-frequency  tidal  and 
gravity  wave  motions  considered  by  Baisley  et  al.  [1983].  The 
occurrence  of  strong  radar  echoes  is  seen  to  be  highly  corre- 
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Fig.  2.  As  in  Figure  i.  for  June  17,  Note  here  the  ' 3-hour  modulation  of  the  velocity  field  and  of  the  lowest  echoing 

levels. 


lated  with  a  particular  phase  of  this  wave  motion.  The  motion 
field  on  June  17  (see  Figure  2)  exhibits  a  dominant  wave  with 
7^3  hours  up  to  ^  1200  AST.  After  that  time,  echoes  are 
sufficiently  sparse  that  no  wave  motion  can  be  discerned.  As  in 
the  case  of  the  ^  7-hour  wave,  the  phase  of  the  ~  3-hour  wave 
motion  appears  to  be  correlated  with  the  occurrence  of  strong 
echoes  at  lower  levels.  Unlike  the  previous  example,  no  low- 
frequency  motions  are  apparent  on  this  day. 

A  coiiipi^nicntary  data  set  is  provided  by  the  rocket  wind 
and  temperature  observations  obtained  with  the  Nike-Hydac 
instrument  package  launched  at  2051  LT  during  salvo  2  on 
June  15.  This  data  set,  while  not  providing  any  information  on 
the  temporal  evolution  of  the  wind  and  temperature  fields, 
does  provide  comprehensive  data  on  the  vertical  structure  and 
the  amplitudes  of  various  scale«  of  motion  a*  a  particular  time 
over  a  height  range  of  55-130  km.  The  zonal  and  meridional 
wind  profiles  and  the  temperature  profile  obtained  with  the 
Nike-Hydac  are  shown  in  Figures  3  and  4. 

The  dominant  vertical  scales  in  the  two  velocity  profiles  are 
in  the  range  ^5  15  km,  with  smaller-amplitude  motions  (with 
A.  $  2  km)  at  lower  levels  and  large-scale  (/.,  ~  30  km)  and 
large-amplitude  motion  apparent  at  upper  levels.  Overall,  it  is 
the  larger  amplitudes  associated  with  the  larger  vertical  scales 
that  appear  to  be  most  significant  here. 


The  temperature  st.  wn  in  Figure  4,  in  contrast  to  the 

velocity  data,  suggest  ‘  otions  with  both  large  and  small 
vertical  scales  contribui.  .ppreciably  to  temperature  fluctu¬ 
ations,  with  the  maximum  gradients  associated  primarily  with 
the  minimum  vertical  scales.  As  will  be  seen  in  section  3,  it  is 
the  temperature  gradient  that  determines  whether  a  particular 
wave  motion  is  stable  or  (convectively)  unstable. 

Thus  the  velocity  and  temperature  data  provide  us  with  two 
diflTerent  views  of  the  same  wave  field.  It  w'ill  be  shown  in 
section  3.  however,  that  these  two  fields  permit  a  much  more 
complete  specification  of  wave  parameters  and  propagation 
characteristics  than  is  possible  with  either  field  alone.  In  some 
instances,  in  fact,  the  two  sets  of  data  are  redundant  and 
constitute  a  test  of  the  experimental  techniques. 

3.  Gravity  Wave  Structure 
The  purposes  of  this  section  are  to  review  linear  gravity 
wave  structure  and  to  demonstrate  what  wave  parameters  can 
be  inferred  using  the  radar  and  rocket  velocity  and  temper¬ 
ature  data  collected  during  the  STATE  experiment.  The  re¬ 
sults  will  then  be  employed  in  section  4  to  examine  the  struc¬ 
ture  and  evolution  of  the  gravity  wave  field. 

Throughout  much  of  the  middle  atmosphere,  the  mean  state 
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Fig.  3.  Zona!  and  meridional  velocity  profiles  derived  from  the 
three-axis  accelerometer  data  obtained  during  Salvo  2  on  June  15. 


and  vertical  wavelengths  which  satisfy 

/..■«(47tH)*  (3) 

where  /  is  the  inertial  frequency,  k[  =  2n  /,)  is  horizontal  wave 
number,  c  is  phase  speed,  li  is  the  mean  flow  in  the  direction  of 
wave  propagation,  and  H  is  the  atmospheric  scale  height. 
Then  the  vertical  wave  number  is  simply 

m  =  2n;/..  =  ±  /V  (c  —  u)  (4) 

Thus  provided  the  mean  atmospheric  stratification  is  known, 
measurement  of  the  vertical  wavelength  of  a  gravity  wave 
implies  knowledge  of  the  intrinsic  phase  speed,  c  -  u,  of  the 
motion. 

Now  assuming  a  wave  motion  of  the  form 

d>'  =  (fig  exp  [/{A:.v  -I-  m~  —  ket)]  (5) 

the  linear  thermodynamic  energy  and  continuity  equations 
yield 

d.'/fl.  =  i<'/(ii  —  c)  (6) 

Thus  a  wave  motion  achieves  a  state  of  convective  instability 
(where  0.'  g  —0.)  at  the  level  of  maximum  horizontal  velocity 
in  the  direction  of  wave  propagation.  For  motions  with  a)  » /, 
this  corresponds  to  the  condition 


is  characterized  by  a  large  Richardson  number. 


u.'  >  N 


(7) 


Rig  = 


(It/  +  r,2) 


»  I 


(I) 


where  N.  u.,  and  v.  are  the  Brunt-Vaisala  frequency  and  the 
mean  zonal  and  meridional  wind  shears,  respectively.  Under 
these  conditions,  gravity  wave  parameters  vary  little  over  a 
vertical  wavelength  and  are  described  well  by  WKB  theory.  In 
addition,  we  assume  intrinsic  frequencies  (i.e..  those  measured 
relative  to  the  local  mean  flow)  such  that 


It  must  be  noted,  however,  that  the  region  of  instability  is 
expected  where  U.'  is  a  minimum  and  it'  is  a  maximum,  rather 
than  in  regions  of  maximum  shear. 

For  gravity  waves  with  intrinsic  frequencies  near  /.  the  Co¬ 
riolis  acceleration  becomes  important,  and  the  expression  for 
vertical  wave  number  becomes 


k^{N^  —  io^)  k^N^ 
i(0^  -f^)  "  (tu^  -f^) 


(8) 


/  ^  «  (u^  =  k\c  —  u)^  «  (2)  Another  consequence  of  rotation  is  a  horizontal  perturbation 

velocity  transverse  to  the  direction  of  wave  propagation  given 
by 


T  (“K) 


Fig.  4.  Temperature  profile  derived  from  the  three-axi.s  accelero¬ 
meter  data  obtained  during  Salvo  2  on  June  15.  Note  the  conspicuous 
~  30-km  wavelength  at  upper  heights. 


if  , 

r  =-u'  (9) 

which  approaches  u'  in  magnitude  as  m— » f.  This  altered  wave 
structure  results  in  enhanced  velocity  shears  and  suppressed 
thermal  fluctuations  for  CO  implying  that  such  motions 
may  become  dynamically  unstable  at  amplitudes  smaller  than 
those  required  for  convective  instability  [see  Friits  and  Ras- 
togi,  1985]. 

Regions  of  wave  field  instability  in  the  middle  atmosphere 
occur  in  response  to  several  factors,  including  amplitude 
growth  due  to  the  decrease  in  atmospheric  density  with  height, 
a  reduction  in  the  intrinsic  phase  speed  due  to  propagation  in 
a  sheared  environment,  and  a  reduction  of  vertical  wavelength 
duo  to  -in  increase  in  atinosphciic  Siability  iBalsley  et  at.. 
1983]  or  as  a  result  of  wave  superposition.  The  resulting  insta¬ 
bilities  produce  turbulence,  which  acts  to  restrain  wave  ampli¬ 
tudes.  This  contributes  to  the  diffusion  of  heat  and  constit¬ 
uents  and  causes  an  acceleration  of  the  local  mean  flow  due  to 
a  divergence  of  the  momentum  flux  of  the  gravity  wave  mo¬ 
tions. 

The  gravity  wave  structure  outlined  above  indicates  that  a 
number  of  important  wave  parameters  can  be  inferred  from 
the  rocket  and  radar  data  collected  during  the  STATE  experi- 
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Fig.  5.  DifferenctJ  temperature  and  velocity  profiles.  Note  that  the  largest  gradients  are  associated  with  small  vertical 

wavelengths. 


merit.  The  vertical  wavelength  can  be  measured  using  either 
rocket  or  radar  data  and  also  provides  an  estimate  of  the 
intrinsic  phase  speed,  c  —  u,  using  (4).  Radar  and  rocket  veloc¬ 
ity  data  can  be  used  to  infer  horizontal  perturbation  velocities, 
directions  of  propagation,  and,  together  with  the  intrinsic 
phase  speed,  wave  amplitudes  relative  to  saturation  values. 
The  radar  data  also  yield  estimates  of  kc  and  u,  and  thus  of  a» 
and  using  (4)  and  (9).  Finally,  the  rocket  temperature  data 
facilitates  estimates  of  O'  and  0^',  permitting  a  check  on  the 
wave  amplitudes  inferred  from  the  velocity  data.  Thus  subject 
to  reasonable  assumptions,  all  of  the  important  parameters 
can  be  inferred  for  gravity  waves  that  are  sufficiently  well 
defined  in  the  STATE  data.  Several  analyses  which  examine 
the  evolution  of  the  gravity  wave  field  and  its  role  in  the 
generation  of  turbulence  are  presented  in  the  following  sec¬ 
tions. 

4.  Analysis  of  Salvo  2  Wave  Field  Data 

The  most  comprehensive  data  on  the  thermal  and  velocity 
fluctuations  due  to  atmospheric  wave  motions  obtained  in  the 
STATE  experiment  were  those  collected  during  salvo  2  on 
June  15.  These  data  include  radar-inferred  horizontal  and  ver¬ 
tical  velocities  between  ~80  and  90  km,  with  temporal  and 
vertical  resolutions  of  ~2  min  and  300  m  (see  Figure  1),  as 
well  as  rocket-inferred  velocity  and  temperature  data  between 
55  and  130  km  with  a  height  resolution  of  ~250  m  (see  Fig¬ 
ures  3  and  4).  As  noted  in  section  2,  these  data  display  a  range 
of  wave  moiions  occurring  on  various  temporal  and  spatial 
scales.  The  analysis  presented  in  this  section  will  focus  on 
those  motions  with  vertical  wavelengths  of  ~2,  12,  and  30  km. 

4.1.  The  /.j  ~  30-km  Wave 

We  assume  here  thr’>  the  fluctuations  noted  in  the  rocket 
and  ladar  data  presented  in  Figures  1-4  are  a  manifestation  of 
internal  gravity  waves,  in  order  to  apply  the  results  of  section 
3.  This  is  almost  certainly  true  for  the  small-  and  medium- 
scale  motions  (/.^  <  20  km).  But  for  larger-scale  motions  (/,  ~ 
30  km),  there  is  a  possibility  that  the  fluctuations  are  due  in 
part  to  tidal  motions,  with  semidiurnal  and  terdiurnal  compo¬ 
nents  contributing  preferentially.  Because  of  the  enormous 
daily  variability  of  the  low-frequency  wave  structure  in  the 
radar  velocity  field  (there  is  no  evidence  of  a  low-frequency 


motion  on  June  17  in  Figure  2)  and  the  general  consistency  of 
the  radar-  and  rocket-inferred  wave  parameters  for  the  large- 
scale  motion  obtained  in  this  section,  however,  we  assume  that 
gravity  waves  are  the  major  contributor  to  the  motion  field 
near  and  above  the  mesopause.  We  recognize,  nevertheless, 
that  inferred  wave  parameters  are  likely  influenced  by  tidal 
motions  at  upper  levels  to  some  degree. 

By  far  the  dominant  wave  structure  apparent  in  the  rocket 
velocity  and  temperature  profiles  is  that  with  A.  ~  30  km,  seen 
at  upper  levels  in  Figures  3  and  4.  In  order  to  obtain  a  reason¬ 
able  fit  to  these  data,  we  first  difference  each  data  set  at  a 
250-m  interval  to  remove  the  mean  and  linear  trend.  It  is 
particularly  important  to  remove  the  trend  from  the  tempier- 
ature  data  because  of  the  rapid  increase  of  temperature  with 
height  above  the  cold  summer  mesopause.  The  differenced 
profiles  are  displayed  in  Figure  5.  As  noted  in  section  2,  the 
maximum  temperature  gradients  are  clearly  associated  with 
the  smallest  vertical  wavelengths. 

We  assume  a  saturated  (constant  amplitude)  wave  motion 
with  uniform  structure  above  85  km.  This  is  appropriate  be¬ 
cause  N  and  (c  -  u)  do  not  change  significantly  in  a  vertical 
wavelength  (Ri'o  »  !)•  We  then  fit  the  differenced  velocity  and 
temperature  data  only  above  the  mesopause  to  insure  that  the 
temperature  minimum  at  that  height  does  not  bias  the  ampli¬ 
tude  or  phase  of  the  inferred  temperature  wave.  The  least 
squares  fit  to  the  diflferenced  temperature  data  occurs  for  = 
30  km  and  is  shown  in  Figure  6.  The  temperature  amplitude 
of  this  wave  is  59.7'C,  with  a  maximum  occurring  at  107.9  km. 

Similar  fits  to  the  differenced  zonal  and  meridional  velocity 
profiles  were  performed  and  the  results  tabulated  in  Table  1. 
The  most  notable  feature  of  the  velocity  field,  also  apparent  in 
Figure  3,  is  that  the  zonal  and  meridional  velocities  are  nearly 
in  quadrature  and  of  comparable  amplitude,  suggesting  a  low- 
frequency  motion  with  downward  phase  progression  (and 
upward  energy  propagation).  It  must  be  cautioned,  however, 
that  there  are  potential  errors  in  the  amplitude  and  phase  of 
the  two  velocity  components,  particularly  the  zonal,  because 
of  possible  variations  of  the  mean  flow  on  the  scale  of  the 
wave  motion. 

With  both  temperature  and  velocity  perturbations  specified, 
we  are  in  a  position  to  estimate  the  direction  of  propagation 
of  the  wave  motion  as  well.  This  can  be  done  either  by  deter¬ 
mining  the  direction  of  the  major  axis  of  the  velocity  ellipse  or 
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Fig.  6.  Differenced  temperature  data  and  best  fit  to  large-scale  (30 
km)  wave  motion.  The  amplitude  of  this  motion  is  less  than  that 
required  for  convective  instability. 

from  the  relative  phase  of  the  temperature  and  velocity  fluctu¬ 
ations.  From  the  heights  of  maximum  amplitude  for  the  two 
velocity  components,  we  infer  a  phase  difference  of  82°.  This 
implies  a  maximum  velocity  of  82.5  ms"'  at  91.5  km,  corre¬ 
sponding  to  a  direction  of  propagation  of  either  6°  or  186°, 

Alternatively,  we  anticipate  that  the  wave  motion  will  have 
a  maximum  negative  temperature  gradient  at  the  height  at 
which  the  wave  velocity  is  in  the  direction  of  propagation. 
This  occurs  at  85.4  and  115.4  km,  where  the  inferred  wave 
velocity  is  at  305°.  Thus  on  the  basis  of  the  rocket  data  alone, 
there  is  considerable  uncertainty  in  the  direction  of  wave 
propagation  because  of  the  errors  inherent  in  the  wave  veloci¬ 
ty  determinations  and  possible  contamination  by  tidal  modes. 
It  is  shown  later,  however,  that  a  much  more  consistent  pic¬ 
ture  of  the  wave  motion  emerges  when  the  radar  data  are 
considered  as  well. 

It  was  shown  previously  that  the  ~  30-km  wave  present 


in  the  rocket  data  corresponds  to  a  low-frequency  motion 
with  a  velocity  field  that  rotates  clockwise  with  height.  This 
same  character  is  observed  in  the  radar  velocity  field  shown  in 
Figure  1.  There  is  also  evidence  that  these  velocities  rotate 
clockwise  in  time,  providing  an  additional  suggestion  of  a 
low-frequency  wave  motion.  A  least  squares  fit  to  the  velocity 
field  between  1200  and  2400  LT  at  each  height  was  obtained 
for  a  wave  period  of  7  hours.  The  phases  and  amplitudes 
obtained  for  the  "east"  and  “north”  beams  at  each  height  are 
shown  in  Figure  7.  These  exhibit  a  clear  phase  progression 
with  height,  with  an  average  phase  difference  of  ~  84°  and  an 
"east”  velocity  approximately  40%  larger  than  the  “north” 
velocity.  There  is  also  a  tendency  for  each  velocity  component 
to  increase  in  magnitude  by  —30%  over  the  —  6-km  height 
interval,  a  factor  much  smaller  than  anticipated  in  the  absence 
of  dissipation  for  a  scale  height  of  —4  km.  Using  these  com¬ 
ponent  velocities,  we  determine  the  velocity  ellipse,  and  thus 
the  wave  amplitude  and  the  apparent  direction  of  propaga¬ 
tion,  at  each  height.  Because  the  mean  velocity  has  been  re¬ 
moved  at  each  level,  this  analysis  yields  a  direction  of  propa¬ 
gation  that  is  quite  uniform  with  height  and  more  reliable 
than  that  obtained  using  the  rocket  velocity  and  temperature 
data.  The  inferred  direction  of  propagation  is  222°  +  23°,  ap¬ 
proximately  midway  between  those  values  determined  using 
the  rocket  data. 

The  radar  configuration  also  provides  information  on  the 
vertical  velocities  due  to  gravity  wave  motions.  For  low- 
frequency  motions,  however,  these  velocities  are  of  the  order 
of  the  system  resolution  used  during  the  STATE  expieriment. 
Because  this  results  in  significant  errors  in  the  estimates  of 
amplitude  and  phase,  these  data  are  not  used  in  this  analysis. 

We  now  proceed  to  calculate  the  ratio  of  observed  wave 
amplitude  relative  to  that  required  for  monochromatic  wave 
saturation  [see  Fritts  and  Rastogi,  1985].  From  (6),  this  ratio 
may  be  expressed  in  terms  of  the  wave  velocity  or  temperature 
as 

u  T' 

a  = - =  ^ — -  (10) 

(m  -  c)  {T^  +  g/c^) 

Above  the  mesopause, 

T; -b  g/c,  «  16°C/km  (11) 


TABLE  1.  Wave  Parameters  Inferred  From  Rocket  and  Radar 
Data  for  Large-Scale.  Low-Frequency  Wave  Motion  Present 
During  Salvo  2  on  June  15.  1983 


Parameter 

Rocket 

Radar 

Data  height  interval,  km 

85-130 

83-89 

A.,  km 

30 

~70 

Period,  hours 

~I7 

T.  A  (°C)/height  of  max,  km 

59.7/107.9 

... 

A  (m  s-’)/height  of  max,  km 

42.2/98.1 

40  ±  4/"- 
(beam  RXl) 

A  (m  s' 'i/height  of  max,  km 

82.3/91.3 

28  ±  3/--- 
(beam  RX2) 

Adi,  deg 

82 

84 

A  (m  s  '(/height  of 
max,  km 

82. 5, '91. 5 

40  +  4/--- 

Direction  of  'uimax-  dc8 

6  or  186 

222  ±  23 

Direction  from  T’  phase,  deg 

305 

Degree  of  saturation  u'lT.' 

0.76/0.78 

0.42 

Rotation  with  z 

clockwise 

clockwise 

Phase  progression 

downward 

downward 

u',  v'  (ms"’) 


Fig.  7.  Phases  and  amplitudes  of  the  horizontal  wind  components 
for  the  7-hour  wave  inferred  from  the  radar  data.  Beams  RXl  and 
RX2  are  oriented  26  north  of  east  and  west  of  north,  respectively. 
Phase  increases  with  time. 
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zontal  Wdvelength  mav  be  estimated,  using  (2)  as 


|T,'|  =  mr  =  12.5  C/km  (12) 

yielding  an  estimate 

a  =  0.78  (13) 

Alternatively,  with  (11)  and  a  mean  temperature  ~300’C,  we 
obtain 

N"  =|:  [t  +  9/Cp]  =  5.2  X  10^'*  (14) 

Together  with 

|u'L..  =  82.5ms-'  (15) 

and  (4),  this  implies 

Ic  —  iJ|  =  109  ms'  (16) 

and 


a  =0.76  (17) 

Thus  despite  the  obvious  uncertainties  in  our  estimates  of  the 

wave  parameters  using  the  rocket  temperature  and  velocity 
profiles,  these  independent  estimates  of  wave  amplitude  are  in 
remarkably  good  agreement. 

A  similar  exercise  using  the  radar  data  available  at  lower 
heights,  with  i  17  km  and  =  1.05  x  10"^  s'^  (due  to 
the  much  smaller  mean  temperature),  yields 

Ic  -  u|  =  88  m  s'*  (18) 

and 

a  =  0.42  (19) 

The  amplitude  of  the  wave  motion  is  smaller  here,  consistent 
with  the  observed  growth  with  height  in  the  radar  data.  Also, 
the  vertical  wavelength  is  compressed  relative  to  that  value 
inferred  at  upper  levels  because  of  the  much  larger  stratifi¬ 
cation  resulting  from  a  small  mean  temperature  and  large 
mean  temperature  gradient  immediately  above  the  mesopause. 
It  was  this  response  to  the  thermal  structure  that  led  Baisley  et 
ill.  [1983]  to  suggest  that  the  summer  mesopause  should  be  a 
preferred  site  for  gravity  wave  and  tidal  saturation  and  turbu¬ 
lence  generation. 

With  the  available  radar  data,  we  are  also  in  a  position  to 
estimate  the  intrinsic  period  and  the  horizontal  wavelength  of 
the  wave  motion.  Using  (9)  and  the  magnitudes  of  the  pertur¬ 
bation  velocities  parallel  and  perpendicular  to  the  direction  of 
wave  propagation,  we  find  that 

2n  2tr|r'| 

T  =  —  =  — ;  a:  8.8  hours  (20) 

w  f\u'\ 

The  intrinsic  period  can  also  be  estimated  from  the  measured 
period,  the  intrinsic  phase  speed  inferred  from  the  vertical 
wave  structure,  and  the  magnitude  and  direction  of  the  mean 
flow,  32  m  s"  *  at  242  /  This  yields  u  a:  30  m  s"  *  and 

Ti  =  Tc/ic  —  li)  a  9,4  hours  (21) 

As  in  the  case  of  wave  amplitude,  the  consistency  of  these 
estimates  suggests  that  the  wave  is  described  well  by  the  radar 
data.  With  the  intrinsic  phase  speed  and  frequency,  the  hori- 


/. 


27t(f  —  u) 


a  3  X  lO*  km 


C) 


(22) 


As  noted  earlier,  this  low-frequency  wave  motion  has  sev¬ 
eral  characteristics  in  common  with  anticipated  tidal  motions, 
including  a  near-tidal  period,  a  fairly  large  vertical  wavelength 
(z.  >  17  km),  and  a  cIockwi.se  rotation  of  the  velocity  field 
with  height  and  lime.  We  believe,  however,  that  the  very 
sudden  appearance  of  the  wave  motion,  together  with  the  con¬ 
sistency  of  the  inferred  wave  structure,  suggests  that  the  wave 
is  not  a  tidal  motion.  Our  determination  of  the  wave  parame¬ 
ters  may.  nevertheless,  have  been  affected  to  some  degree  by 
the  presence  of  tidal  components.  We  also  note  that  the  effects 
of  this  wave  motion  appear  to  be  virtually  identical  to  those 
attributed  to  tidal  and  low-frequency  gravity  wave  motions 
near  the  summer  mesopause  by  Baisley  et  al.  [1983]. 

We  now  examine  the  consequences  of  a  saturating  wave 
field  dominated  by  this  30-km  wave  motion.  From  the  analy¬ 
sis  presented  by  Fritls  and  Rastogi  [1985],  we  see  that  a  wave 
with/'fu  =  0.72  is  dynamically  unstable  for 

a  g  0.82  (23) 

Thus  the  30-km  wave  al  upper  levels  is  nearly  dynamically 
unstable,  even  in  the  absence  of  other  wave  motions.  There  is 
little  doubt  that  the  large-scale  (z.  ^  5  km)  wave  field  is  dy¬ 
namically  and/or  convectively  unstable  at  some  locations 
when  all  components  are  considered. 

At  the  level  of  radar  observations,  the  7-hour  wave  has  an 
amplitude  substantially  less  than  that  required  for  instability. 
This  does  not  imply  that  the  overall  wave  field  is  stable,  how¬ 
ever.  Indeed,  wave  amplitudes  smaller  than  monochromatic 
saturation  values  have  been  inferred  in  several  observational 
studies  of  gravity  wave  saturation  [Smirfi  and  Fritts,  1983; 
Meek  et  al..  1985a;  Reid  and  Vincent,  1987a],  and  the  numeri¬ 
cal  study  by  Fritts  [1985]  indicates  that  saturation  may  occur 
at  small  wave  amplitudes  due  to  a  superposition  of  waves. 

Because  the  velocity  field  in  the  last  -- 12  hours  of  Figure  1 
clearly  is  dominated  by  the  7-hour  wave,  it  seems  appropriate 
to  ask  whether  this  motion  is  related  to  the  occurrence  of 
strong  radar  echoes  and  to  the  turbulence  enhancements 
thought  to  be  responsible  for  these  echoes.  To  answer  this 
question,  we  have  shown  contours  of  S/N  above  0  dB  for  this 
period  in  Figure  8.  As  suggested  by  the  velocity  field  in  Figure 
I.  regions  of  maximum  S/N  appear  as  descending  layers  about 
7  hours  apart,  indicating  a  strong  correlation  with  the  7-hour 
wave  motion.  This  correlation  is  further  enhanced  if  we  plot 
the  phase  of  the  wave  motion  corresponding  to  the  most  un¬ 
stable  portion  of  the  wave  field  anticipated  from  the  dis¬ 
cussion  in  section  2.  These  are  shown  by  the  solid  lines  in 
Figure  8  and  indicate  the  locations  in  the  wave  field  where  the 
temperature  gradient  is  most  negative  and  where  the  pertur¬ 
bation  velocity  is  in  the  direction  of  wave  propagation  (222°), 
as  seen  in  Figure  1.  Thus  the  radar  data  provide  strong  evi¬ 
dence  of  a  significant  enhancement  in  turbulence  activity  at 
precisely  that  phase  of  the  wave  motion  where  instability  and 
dissipation  are  expected  to  occur  on  theoretical  grounds.  Fur¬ 
ther  evidence  of  this  correlation  is  provided  by  the  findings  of 
Watkins  et  al.  [1987],  in  which  the  height  of  the  maximum 
inferred  energy  dissipation  rate  was  observed  to  descend  with 
time  and  to  coincide  with  the  maximum  in  S/N  in  Figure  8. 
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Fig.  8.  Time-height  cross  seciion  of  radar  S  N  for  June  IS,  from 
1000  to  2400  LT.  The  slanted  lines  denote  the  most  unstable  phase  of 
the  7-hour  wave  inferred  from  the  temperature  and  velocity  data  and 
correlate  well  with  regions  of  maximum  S/N. 

This  observation  has  significant  implications  for  the  diffusion 
of  heat  and  constituents  and  suggests  a  large  effective  Prandtl 
number  at  these  heights  [C/wo  and  Schoeberl,  1984;  Fritts  and 
Dunkenoti.  1985;  Strobe!  et  al..  1985].  ft  should  again  be 
stressed  that  this  phase  of  the  wave  motion  corresponds  to  a 
minimum  of  the  perturbation  velocity  shear. 

Returning  to  the  differenced  rocket  temperature  and  veloci¬ 
ty  data  (Figures  5  and  6).  we  see  that  there  is  a  tendency  for 
enhanced  small-scale  wave  structure  (/..  ~  2-3  km)  to  occur 
near  those  heights  where  the  temperature  gradient  associated 
with  the  7-hour  wave  is  most  negative  ( ^  85  and  1 15  km).  This 
suggests  that  large-scale  wave  motions  act  as  a  source  of  both 
turbulence  and  small-scale  wave  activity  preferentially  in  re¬ 
gions  where  the  wave  field  is  most  unstable.  In  the  present 
case,  it  appears  that  this  is  a  result  of  a  dynamical  instability 
of  the  flow,  as  described  earlier  in  this  section.  This  conclusion 
is  supported  by  the  lack  of  large-scale  superadiabatic  lapse 
rates  in  the  rocket  temperature  data  (see  Figure  6).  In  other 
cases  of  gravity  wave  saturation  involving  wave  motions  with 
high  intrinsic  frequencies,  it  appears  that  turbulence  and  wave 
excitation  occur  primarily  as  the  result  of  a  convective  insta¬ 
bility  [Fritts,  1984a;  Fritts  and  Rastogi.  1985].  These  observa¬ 
tions  are  in  qualitative  agreement  with  numerical  studies  sug¬ 
gesting  that  both  dynamical  and  convective  instabilities  may 
lead  to  the  excitation  of  other  internal  gravity  wave  motions 
[Fritts,  1984b,  1985;  Chimonas  and  Grant,  1984].  Finally,  an¬ 
other  process  that  may  be  important  in  the  evolution  of  a 
saturating  gravity  wave  spectrum  [fVeinstock,  1982],  but 
which  appears  not  to  prevent  the  occurrence  of  wave  field 
instability  [Fritts,  1985;  Smith  el  al.,  1987],  is  the  nonlinear 
interaction  of  waves  that  have  attained  large  amplitudes. 

4.2.  The  —  2-  and  I2-km  Waves 

As  is  apparent  from  the  discussion  in  section  2  and  Figures 
1-4,  the  rocket  and  radar  data  collected  during  salvo  2  con¬ 
tain  evidence  of  a  wide  range  of  wave  motions  occurring  at 


various  temporal  and  spatial  scales.  1  he  focus  earlier  in  this 
section  was  on  the  30-km  wave  present  at  upper  levels  becau,se 
of  the  domiiiu.it  role  it  obviously  played  in  the  dynamics. 
However,  other  motions  are  likely  to  be  equally  important  at 
lower  levels  where  the  30-km  wave  has  not  approached  a 
saturated  amplitude. 

One  such  motion,  occurring  at  the  height  of  the  radar  echo 
and  appearing  in  both  zonal  and  meridional  rocket  velocity 
profiles  (see  Figure  3),  is  found  to  have  a  least  squares  vertical 
wavelength  of  A,  =  12  km  in  the  interval  from  75  to  95  km.  As 
in  the  case  of  the  30-km  wave,  this  motion  is  assumed  to  be 
saturated  and  thus  to  have  a  constant  amplitude  with  height. 
The  parameters  estimated  for  this  12-km  wave  are  listed  in 
Table  2. 

The  zonal  and  meridional  velocity  components  are  seen  to 
be  nearly  in  phase,  suggesting  a  high-frequency  motion  of 
horizontal  magnitude  19.9  ms'  propagating  toward  41'.  A 
similar  fit  to  the  temperature  data  yielded  an  amplitude  of 
7.6°C  and  a  phase  which,  together  with  the  velocity  data, 
suggested  a  propagation  toward  358’.  It  must  be  noted,  how¬ 
ever,  that  the  estimate  of  temperature  amplitude  is  subject  to 
considerable  error  because  of  the  mean  thermal  structure  near 
the  mesopause.  Nevertheless,  we  estimate  the  wave  amplitude 
relative  to  the  saturated  value  as  before,  assuming  =  6.4 
X  10  s“^,  and  obtain  ratios  of  0.41  and  0.77  from  the 

velocity  and  temperature  data,  respectively.  In  this  case,  the 
stratification  experienced  by  the  high-frequency  motion  is  as¬ 
sumed  to  be  that  due  to  the  mean  state  and  the  7-hour  wave 
motion.  These  results  suggest  that  this  motion  may  contribute 
significantly  to  the  attainment  of  a  saturated  wave  field  at  the 
heights  of  the  radar  echoes.  In  particular,  this  high-frequency 
motion  may  be  the  cause  of  some  of  the  rapidly  descending 
S/N  maxima  embedded  within  the  dominant  7-hour  periodic¬ 
ity  seen  in  Figure  8.  These  maxima  suggest  a  wave  period  of 
~  1  hour,  but  because  of  the  large-amplitude  7-hour  motion, 
no  such  period  was  conspicuous  in  the  radar  data. 

One  final  wave  motion  in  the  salvo  2  data  that  warrants 
some  discussion  is  the  very  conspicuous  ~  2-km  structure 
seen  in  the  temperature  gradient  data  in  Figure  6.  The  small 
vertical  wavelengths  (and  small  intrinsic  phase  speeds)  of  these 
motions  suggest  that  they  must  have  been  generated  locally. 
Additionally,  the  very  large  (convectively  unstable)  amplitudes 
occurring  near  the  minima  of  the  large-scale  dT/dz  profile 
imply  that  such  motions  are  generated  preferentially  where  the 


TABLE  2.  Wave  Parameters  Inferred  From  Rocket  Data  for  the 
A.  ~I2  km  Wave  Motion  Observed  at  Intermediate  Heights 
During  Salvo  2  on  June  15,  1983 


Parameter 

Rocket 

Data  height  interval,  km 

75-95 

A.,  km 

12 

Period  (intrinsic) 

short 

T' .  A  (°C)/height  of  max,  km 

ii'.  A  (m  s  'i/height  of  max,  km 

13.4/89.4 

i  ',  A  (m  s  'l/height  of  ma;;.  km 

15.2/88.7 

AJi,  deg 

22 

Ih'1.  a  (m  s  '(/height  of  max,  km 

19.9/89.0 

Direction  of  deg 

41 

Direction  from  T'  phase,  deg 

358 

Degree  of  saturation  u'lT.' 

0.41/0.77 

Wave  motion  was  not  distinguishable  from  the  low-frequency 
motion  in  the  radar  data. 
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large-scale  wave  field  is  unstable.  Because ‘the  small-scale  mo¬ 
tions  are  themselves  forced  to  unstable  amplitudes,  they  con¬ 
tribute  to  the  generation  of  turbulence  and  the  cascade  of 
wave  energy  to  smaller  scales 

5  Anai  vsisof  Jim  17  Radar  Data 

As  described  in  section  2,  the  radar  velocity  field  obtained 
on  June  P  exhibits  a  pronounced  wave  motion  and  a  corre¬ 
sponding  modulation  of  the  echo  occurrence  (and  intensity)  at 
lower  levels  w  ith  a  period  of  ^  3  hours  (see  Figure  2).  This  is 
very  similar  to  (he  correlation  between  the  velocity  field  and 
echo  occurrence  observed  for  the  7-hour  wave  during  salvo  2 
and  analyzed  in  section  4.  The  purpose  of  this  section  is  to 
determine  whether  the  enhanced  echoes  can  again  be  related 
to  the  most  unstable  portion  of  the  wave  field. 

A  least  squares  fit  to  the  velocity  data  between  ^  83  and  88 
km  from  0400  to  1000  LT  was  obtained  for  a  wave  period 
T  190  min.  This  wave  motion  is  also  observed  above  and 
below  this  height  interval,  but  a  lack  of  sufficient  echo  leads  to 
a  reduction  in  the  number  of  accepted  data  points  and  a 
degradation  of  the  fit  at  these  heights.  The  phase  of  the 
motion  in  the  two  radar  beams  is  shown  as  a  function  of 
height  in  Pigurc  9.  From  these  data  we  infer  a  vertical  wave¬ 
length  of  ^  38  km  and  a  corresponding  intrinsic  phase  speed 
of  ~  150  m  s'‘,  though  there  is  some  uncertainty  in  these 
estimates  due  to  the  dilTerent  phase  slopes  in  the  two  beams. 
The  magnitude  of  the  velocity  fluctuation  and  the  inferred 
direction  of  propagation  are  show  n  as  functions  of  height  in 
Figure  10.  These  data  suggest  an  increase  of  wave  amplitude 
with  height  from  ^  20  m  s  “  ‘  at  83  km  to  ~  26  m  s  “ '  at  86  km 
and  a  direction  of  propagation  of  ~  205  .  The  increase  in 
wave  amplitude  with  height  is  less  than  that  required  for  con¬ 
servative  wave  motion,  suggesting  that  the  wave  is  likely  to  be 
dissipating  and  contributing  to  the  occurrence  of  the  radar 
echoes  at  these  heights. 

The  mean  flow  inferred  between  84  and  86  km  is  ^  60  m 
s  '  at  246  ,  The  component  of  this  flow  in  the  direction  of 
wave  propagation  is  u  ^  45  m  s' which,  with  the  estimated 
intrinsic  phase  speed,  implies  a  phase  speed  c  -  195  m  s”' 
and  a  horizontal  wavelength 

=  cT  -  2230  km  (24) 


PHASE  (rad)  u',  v’  (ms"’) 


Fig.  9.  Phases  and  amplitudes  of  the  horizontal  wind  components 
for  the  190-min  wave  seen  in  the  radar  data  on  June  17.  Phase  in¬ 
creases  with  time. 


AMPLITUDE  (ms-'J  DIRECTION  (°) 


Fig.  10.  Magnitude  and  direction  of  the  horizontal  motion  for  the 
190-min  wave.  Note  that  the  direction  of  propagation  corresponds  to 
the  perturbation  velocity  at  levels  where  the  S  N  is  enhanced 

This  suggests  an  intrinsic  period 
27r 

T  = - ~  247  min  (25) 

ktc  -  li) 

It  must  be  recognized,  however,  that  these  values  reflect  the 
uncertainty  in  our  estimate  of  the  vertical  wavelength  of  the 
wave  motion  and  are  likely  lO  be  less  accurate  than  those 
values  estimated  using  rocket  and  radar  data  in  section  4. 

Thus  we  have  again  obtained  a  reasonably  consistent  de¬ 
scription  of  the  wave  field  on  the  basis  of  the  available  radar 
velocity  data.  This  view  is  reinforced  by  noting  that  the  phase 
of  the  wave  motion  in  which  the  S  N  (and  presumably  insta¬ 
bility  and  turbulence)  is  enhanced  corresponds,  as  in  the  case 
of  the  7-hour  wave,  to  that  whicu  is  expected  to  be  most 
unstable  on  the  basis  of  the  discussion  in  section  3.  This  re¬ 
sults  in  perturbation  velocities  at  the  lowest  echoing  levels 
that  are  in  the  direction  of  wave  propagation,  suggesting  an 
easy  way  to  infer  this  quantity  in  other  data. 

6.  Summary  and  Conc  lusions 

We  have  pre.sented  analyses  of  several  gravity  wave  motions 
observed  in  the  mesosphere  and  lower  thermosphere  using 
data  collected  during  the  STATE  experiment.  This  expieriment 
was  designed  to  examine  the  conditions  responsible  for  the 
strong  Poker  Flat  MST  radar  echoes  occurring  near  the 
summer  mesopause.  The  rocket  and  radar  data  obtained  have 
provided  a  unique  picture  of  both  the  spatial  and  temporal 
structure  of  the  gravity  wave  field  present  at  the  time  of  the 
experiment. 

The  combined  rocket  and  radar  temperature  and  velocity 
data  collected  during  salvo  2  on  June  15  were  found  to  permit 
a  more  complete  specification  of  the  wave  field  than  could 
have  been  accomplished  using  either  system  alone.  Clearly  the 
dominant  motion  was  the  large-scale,  low-frequency  wave  ob¬ 
served  by  both  systems.  Using  the  rocket  data  at  upper  levels 
and  the  radar  data  near  the  mesopause  and  assuming  a  grav¬ 
ity  wave  rather  than  a  tidal  structure,  it  was  possible  to  infer 
all  of  the  principal  wave  parameters  including  amplitude,  hori¬ 
zontal  and  vertical  wavelength,  phase  speed,  intrinsic  period, 
and  direction  of  propagation.  Specifically,  this  wave  motion 
was  found  to  have  an  amplitude  comparable  to  that  required 
for  dynamical  instability  at  upper  levels.  Near  the  mesopause 
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lhai  phase  ol  llic  \sase  that  was  c\pcctcd  ts)  be  most  unstable, 
bases!  on  the  inrerresl  wa\e  struetiire.  was  I'siund  to  sMrrclate 
)er\  well  .Mth  enhaneed  S  N  in  the  radar  data  (and  t>rcsum- 
abl>  with  inereases'  instabihts  and  lurbuleneel.  suggesting  that 
ihe  generation  of  intense  turbulenee  is  highly  le>eali/cd  within 
the  w.ne  lield  and  thus  that  the  Prandtl  number  due  to  wave- 
generated  turbulenee  is  large  [Fniis  uiul  Diinkcrion.  1985], 
This  w,ne  was  estimated  to  propagate  toward  the  southwest. 

\  'eeond  wave  mo’K'ii  with  a  vertieal  wavelength  '  12 
km  was  also  dediieed  using  the  roeket  data  near  the  meso- 
pause  The  amplitude  oT  this  motion  was  estimated  to  be  suf- 
lieient  lo  eontribute  appreeiably  to  the  saturation  of  the  low- 
Irequenev  wave  motion  in  this  regu'ti  .A  high  intrinsic  fre- 
quenev  was  inferred  both  from  the  phase  consistency  between 
the  zonal  .ind  meridional  velocities  and  from  the  high- 
Ireqiiencv  i  /  '  1  hourl  modulation  of  the  radar  S  N  data. 
Tmallv.  this  wave  im'tu'ii  was  found  to  propagate  toward  the 
north 

Also  conspicuous  in  the  rocket  temperature  gradient  data. 
Ill  particular,  were  wave  motions  with  / .  ^  2  km.  These  were 
ob^erved  with  extremelv  large  (tinstablel  amplitudes  near  rc 
glon^  in  which  the  large-scale  wave  field  was  itself  most  unsta¬ 
ble  The  very  slow  vertieal  propagation  and  large  rate  of  dissi¬ 
pation  anticipated  for  such  motums  further  suggest  that  they 
may  result  fri'm  the  instability  (if  the  large-scale  wave  field. 

AiK'ther  low-frequency  event  observed  in  the  radar  data  on 
Itine  1'  was  analyzed  because  it  appeared  to  be  related  to  a 
niiidtilation  of  the  lowest  heights  from  which  useful  radar 
echoes  were  ('btained  This  wave  had  a  period  T  ^  190  min. 
with  clearly  defined  amplitude  and  phase  near  the  mesopause. 
permuting  a  fairly  unambiguous  specification  of  the  wave 
structure.  As  in  the  ease  of  the  7-hour  wave,  a  strong  corre¬ 
lation  between  the  most  unstable  phase  of  the  wave  motion 
■ind  the  occurrence  of  large  S  N  was  observed,  providing  fur¬ 
ther  evidence  that  turbulence  generation  is  localized  within 
the  wave  field.  This  wave  motion  was  found  to  propagate 
toward  the  south-southwest 

The  data  tinalyses  presented  here  suggest  that  occurrences 
of  gr.ivity  wave  saturation  in  the  mesosphere  and  lower  ther¬ 
mosphere  are  largely  consistent  with  theoretical  expectations 
of  saturation  via  a  dynamical  or  convective  instability  of  the 
w.ive  field  [//la/i/es.  1967;  Limizen.  19X1;  Iritl'i  and  Raslinii. 
I9S5]  For  large-scale  motions  there  was  found  to  be  a  strong 
correlation  between  large  radar  S  N  (and  large  inferred  turbu¬ 
lence  intensities)  and  the  most  unstable  phase  of  the  wave 
motion  This  suggests  that  turbulence  generation  and  effects 
may  be  strongly  localized  within  the  wave  field,  perhaps  re¬ 
sulting  in  a  large  efTective  Prandtl  number,  as  anticipated  by 
Chao  and  SLhoeberl  [1984']  and  frills  and  Dunkerlon  [198.“)]. 
Evidence  was  also  found  of  the  excitation  of  other  gravity 
wave  motions  at  small  scales  preferentially  near  regions  of 
maximum  instability  of  the  large-scale  wave  field.  Whether 
these  motions  were  the  results  of  wave-wave  interactions  of 
the  sort  considered  by  H'einslm  k  [1982]  and  Frills  [1985]  or 
a  consequence  of  dynamical  or  convective  instabilities  [Fri/fs. 
19846.  1985]  could  not  be  determined. 

It  should  also  be  noted  that  there  is  a  tendency  to  select  for 
analysis  those  wave  motions  with  the  largest  amplitudes  and 
hence  the  largest  vertical  wavelengths  and  intrinsic  phase 
speeds  These  motions  are  expected,  because  of  their  large 
energies  and  large  vertical  group  velocities,  to  be  the  most 
significant  in  terms  of  turbulence  production.  However,  mo¬ 


tions  with  small  intrinsic  frequencies  are  not  expected  to  con¬ 
tribute  appreciably  to  the  transport  of  momentum  because  of 
their  relatively  small  vertical  perturbation  velocities  [T/iuenf 
and  Reid.  1983;  Frills.  1984a].  Indeed,  those  motions  analyzed 
in  this  paper  were  propagating  toward  the  north  or  the  south¬ 
west  and  thus  could  not  have  contributed  appreciably  to  the 
deceleration  of  the  mean  easterly  flow  in  the  mesosphere.  Be¬ 
cause  of  this  and  other  related  data  [Smiih  and  Fritis.  198.3; 
.Miek  el  at..  19856;  I  ineeni  and  Frills.  1987].  it  would  appear 
essential  to  consider  those  factors  aflecting  the  meridional  as 
well  as  the  zonal  propagation  of  gravity  waves  in  any  compre¬ 
hensive  treatment  of  grav  ity  wave  elTects  in  the  middle  atmo¬ 
sphere.  This  element  of  the  problem  has  been  largely  neglected 
in  most  studies  to  date. 
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